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Abstract 

The thermodynamical properties of a quantized electromagnetic field inside a box with perfectly 
conducting walls are studied using a regularization scheme that permits to obtain finite expressions 
for the thermodynamic potentials. The source of ultraviolet divergences is directly isolated in 
the expression for the density of modes, and the logarithmic infrared divergences are regularized 
imposing the uniqueness of vacuum and, consequently, the vanishing of the entropy in the limit 
of zero temperature. We thus obtain corrections to the Casimir energy and pressures, and to the 
specific heat that are due to temperature effects; these results suggest effects that could be tested 
experimentally. 

PACS numbers: 03.70.+k, ll.lO.Wx, 42.50.Lc 
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I. INTRODUCTION 



The theory of blackbody radiation in cavities has played a decisive role in the development 
of quantum physics. In particular, the existence of a zero-point energy Q with measurable 
physical effects, as predicted by Casimir in 1948 0], followed from its consistency. It is now 
well established that the fluctuations of the quantum vacuum induce an attractive force 
of magnitude (per unit area) 7r^/ic/240L^ between two perfectly-conducting parallel plates 
separated a distance L. The physical importance of the Casimir effect prompted an extensive 



investigation of the fundamental properties of quantum and thermal 



iuctuations for a large 



Nevertheless, a 




variety of systems with particular dielectric and geometric features 
complete self-consistent theory of quantum fluctuations in closed cavities is far from being 
achieved due to the lack of a clear and well established regularization scheme. 

In the simplest case of an ideal conducting parallel-plate configuration at finite temper- 
ature, relevant thermodynamical quantities such as entropy, internal and free energies_^aiid 
pressure forces have been calculated by many authors following different lines jj, \~~ 
In a pioneering work based on general symmetry considerations. Brown and Maclay 
alyzed the structure of the electromagnetic stress-energy tensor T^^^ for two parallel plates 
immersed in a thermal bath at a temperature T. They found that this tensor can be written 
as a sum of three terms: 

T^'^ = T^;'^ + T^^^ + T^^^^ (1) 

where the first term represents the zero-point Casimir stress tensor for a plate separation 
L, the second is due to the black body radiation between widely separated plates, and the 
third term is a correction that vanishes in the limits T — > and L —>■ oo. The resulting 
energy density coincides with the standard Stefan-Boltzmann expression, E/V = T^^^^ = 
[n'^ /15h'^c^){kBT)'^, but there is also a contribution due to boundary effects in the limit 
T — i> oo: the pressure at the plates deviates from the standard blackbody pressure by a 
ierm that increases linearly with the temperature. This additional term is independent of h 
^; in fact, Boyer also derived it within a purely classical formalism, and Mehra Q], and 
Feinberg et al. jll obtained identical results. This linear term has also been found for other 



geometries and materials 



IG, 



There are some controversial issues regarding the Casimir forces between dielectric parallel 



plates. For instance, the Lifshitz formula 



12| that follows from a direct application of 
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Matsubara's frequencies formula is ambiguous in the infrared regime. The fundamental 
issue seems to be whether the transverse electric (TE) component of the electromagnetic 
field contributes to the Casimir force in that regime. The theoretical answer depends on the 
model used to characterize the dielectric response of the material: at low frequencies, the 
TE mode vanishes if the Drude model for the dielectric function is used but this is not 
the case of the plasma model jlJ]. According to the first approach, there must be a large 
thermal correction to the Casimir force between plates with a separation of one micrometer 
o. la.^ tat no such effect Kas been obsc„ed in cxpC.ents QQ. Mo.cove. Been. 



et al. 



17l | have shown that using the Drude model, the entropy turns out to be finite at zero 



temperature and depends on the parameters of the system, thus violating the third law of 
thermodynamics; these authors proposed to solve the problem with the plasma model and 
the Leontovich surface impedance approach, but the app licability of this scheme has been 



3e mipedance ap; 

Q,Q,Q- stii 
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2l| have concluded that some 



questioned by others [18|, ll9|, |20| . Still other authors 
of the discrepancies are due to the fact that, in the low frequency regime, the atomic nature 
of the material invalidates the use of dielectric functions that characterize bulk properties. 
The problem becomes considerably more cornplicated beyond the parallel plates configu- 



ration. In their pioneering work. Case and Chiu [2^ used the fluctuation-dissipation theorem 
to calculate the mean internal energy inside a perfectly conducting cavity filled with a di- 
electric dispersive medium at finite temperature. For an empty cubic cavity, they obtained 
an expression for the internal energy per unit volume similar to Eq. , and they also found 
a deviation from the Stefan-Boltzmann law at low temperatures. However, Case and Chiu 

did not calculate other thermodynamic potentials that also characterize the system. 

n 

In a later work, Ambjorn and Wolfram [Q] studied the general problem of the constrained 
vacuum fluctuations of scalar fields in D-dimensional Euclidean spaces, subject to either 
Dirichlet or Neumann boundary conditions along P orthogonal directions; they calculated 
the partition function of such a configuration with a dimensional regularization technique. 
However, the resulting partition function for a completely closed cavity [D = P) turned 
out to be divergent, as well as thermodynamic quantities such as the free energy and the 
entropy. Ambjorn and Wolfram claimed that these divergences could be removed introducing 
two further regularization schemes, valid at low and high temperatures separately, but they 
gave no explicit expressions for the finite thermodynamic potentials that may be valid at all 
temperatures. A similar problem also arises when dealing with electromagnetic fluctuations 
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23| who studied the asymptotic 



in rectangular closed cavities, as shown by Santos and Tort 
behavior of the free energy at low and high temperatures using a duality inversion formula. 
Such a duality ^mmetry was formerly discovered by Brown and Maclay for the parallel plate 
configuration 5] and studied by Ravndal and Tollefsen 2^ . More recent works along these 
lines are due to Inui 2^ and Cheng [2^. However, none of these authors report the free 
energy of the EM field in a closed and finite form that could be valid for all temperatures and 
box geometries. As a consequence, no quantitative predictions for temperature corrections 
to the Casimir pressures can be deduced from these works. 

The aim of the present paper is to study the thermal effects of a quantum electromagnetic 
field in closed rectangular cavities with perfectly conducting walls. For this purpose, we pro- 
pose a regularization scheme that permits to calculate the thermodynamical variables and 
obtain finite results that can be eventually tested in laboratories. The analogous pro blem 
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in the zero-temperature case has been addressed in several previous works 
the main results being that the Casimir energy and pressures on the cavity plates depend 
strongly on the geometry (they can even change sign according to the relative ratios of the 
faces). In general, the Casimir pressures produce instabilities and, in the particular case of 
a cube, the forces acting on its faces are repulsive In the present paper we generalize 
the previous results to include a thermal bath in the configuration. In Section 2, the regu- 
larization is achieved using a particular representation of the EM mode density that permits 
to isolate the contributions of ultraviolet divergences (arising from zero-point fiuctuations, 
to the internal and free energies). For this purpose, we first isolate the contribution to the 
mode density in a large cavity, as given by Weyl's asymptotic expression, including correc- 
tion terms proportional to the edges of the configuration. Next, the infrared divergences are 
regularized by imposing the uniqueness of the vacuum state, i.e. the state without photons 
of any frequency, in such a way that the entropy of the system at zero temperature is zero by 
construction. The regularized free energy is then used to calculate the pressures, which are 
physical observable quantities in experiments. In Section. 3, our results are compared with 
other approaches in order to elucidate the origin of the difficulties arising in the problem of 
quantum and thermal fiuctuations in rectangular cavities. 
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II. THERMODYNAMIC VARIABLES 



For any given field at thermal equilibrium with its surroundings, the contribution to the 
free energy of a mode of frequency is 

= T ln(l - e-"--/^) + ^nu;^ , (2) 

in standard notation, where the second term is the contribution of the zero-point field (from 
now on, we set ks = 1). The total free energy is a sum over the set of modes: 

p{uk)F{uJk.)duJk , (3) 

where p^uj^) is the density of modes. 

If we consider a rectangular cavity with sizes {k = 1, 2, 3) made of an idealized perfect 
conducting material, the boundary conditions imply the discretization of the frequencies in 
the box: 




The density of modes has the form [22, l2r 

n 

where the prime in the summation indicates that the term with all three indices = 
is to be excluded. This particular form of p{u!) reflects the fact that the normal (parallel) 
components of the magnetic (electric) field vanish on the walls of the cavity, and that 
the electromagnetic field admits two polarization states represented by transverse electric 
(TE) and transverse magnetic (TM) modes. These modes may be derived from massless 
scalar Hertz potentials that satisfy mixed combinations of Dirichlet or Neumann boundary 
conditions on the cavity plates [2^ 



As shown in Ref. j2S|, ultraviolet divergences of the zero-point energy can be isolated. For 
this purpose, the spectral density is written in terms of a summation involving the variable 



2 

C 



1 1/2 

(nioi)^ + (^202)^ + (^303)^ , (6) 



which is conjugate to the frequency variable uj. As shown in Ref. j28|, it is convenient to 
separate the density of modes into two parts: 

p{u;)=p^{u;) + Ap{u;), (7) 



where the first term is 



P°°M = ^^'-^(«l+«2 + 03), (8) 



and the second term can be written in the form 

V -^-^ OOSmiuJ Un) CLi 



^P(-) = ^sl^ ^ "«oo) 

n=— OO n=— OO 

OO ' OO ' 

^2 , ^ as 



COs(t^ Mo„o) - 7r~ "OOn) • (9) 

n=— OO 



The above relations are vahd for all values of u; 7^ 0. 

The variable cwn in the above formula can be interpreted as the position of the image 
charges that generate the field correlation functions with the appropriate boundary condi- 



tions on the conducting walls 



28|,|30|. Moreover, since these variables are conjugate to the 



frequency, lower values of the integer rij's in u-a, correspond to higher values of the frequency, 
and vice versa. Notice also that Eq. (jHl) is not well defined for = because the terms 
related to edges diverge, while in the original expression for p{uj\ Eq. (jSl), the term with 
u; = is absent. 

The spectral density p'^iuj) prevails in the limit Oj — ^ 00, but it is not identical to 
the density of modes in free space because the edge terms are important even in this limit. 
However, just for simplicity, we shall call it "free" or "blackbody" density . The use of Ap{uj) 
instead of the full density of states p{uj) for evaluating an extensive variable is equivalent 
to calculating the difference between the values of that variable in bounded and in "free" 
space. 

The contribution of the zero-point field to either the free or the internal energy associated 
to p°°{uj) is always divergent. However, the temperature dependent part of the free energy 
that corresponds to the blackbody radiation is perfectly finite and is given by 



POO 

F^B^T p~Hlog(l-e- 
Jo 



huj/T 



)duj 



The first term, proportional to the volume of the box, is the usual text-book formula; the 
second term is the contribution of the cavity edges to the free energy jjs^]. It is easy to 
see from a dimensional analysis that any contribution to the free energy proportional to 
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the surface of the walls ajOj should be proportional to T^, and similarly a contribution 
related to the vertices should be proportional to T and independent of the parameters Oj. 
As mentioned above, contributions to the spectral density from faces and vertices of the 



cavity do not appear due to a cancellation of the corresponding TE and TM modes 
In any case, a correction of the Stefan- Boltzmann law, possibly as the one given by Eq. (|lUj). 
has been observed during the development of masers technology js^. 

The contribution of the spectral density Ap{uj) to the zero-point free or internal energy 
is the Casimir term 



AEq = — uj/S.p{ijj)duj 
2 Jo 



n 

n "nOO ^^'^ „ "OnO ^^'^ „ "OOn 

The thermal part of the free energy calculated with Ap{uj) turns out to be free of ultraviolet 
divergences. However, some of the integrals involved do not have a unique interpretation as 
shown in the Appendix. In fact, 

2 

n 

+ ^^^X^ [ai/£;(^^noo) + a2fE{vono) + asfEivoon)] , (12) 

n 

where we have defined 

fvA'") = -(giv) - Kv,e{v] 

V V 

hv,E{v) = h{v) + —Kv^e{v) , 
giy) = coth(t') , 

h{v) = -^(-9{v)) (13) 

V dv\v J 

as functions of the dimensionless variable 

ttT 
a 

(at room temperature /c^T ~ 2.6 x 10^^ eV, and this implies that hc/ixkBT ~ 3/im ). Ky{v) 
and Ke{v) are trichotomic functions of v that can take only the values ±1 and 0, depending 



on the contour of integration selected to evaluate the integrals appearing in the volume and 
edge parts, respectively, of the free energy AF. Notice that in the limit T — > 0, the free 
energy AF is equal to the internal energy AEq independently of the value of K. 

In general, given an integral expression for a physical variable, the selection of an integra- 
tion path corresponds to a particular boundary condition for that variable. In the following, 
we shall impose the limiting condition 

dF 

in order to select an appropriate integration path for a given value of Vn- The entropy S of 
the system being 



S = ~ 



dF 
df 



(15) 



Eq. ()14|) implies the uniqueness of the EM vacuum. The blackbody contribution to the 
entropy is 



coo 



47r^ TT^ 



(16) 



and using the integrals given in the Appendix, the remaining part of the entropy can be 
written as 



AS 



rii 

^ [^fE{Vn00) 



-Kv{Vn) +4:hv{Vn) 



+ vLnhE(VnOo) 



naiT 
Ahc 

n 
n 

-^Yl [^/^(^OOn) + vlo^hEiVoOn) 



DnOj 



(17) 



In order to evaluate the limit T ^ of the above formulas, the summation over discrete 
indices can be replaced by an integration over a continuous dummy variable v. 



S 



1 r°° 



dv 



v^[vh\v) - —,Kv{v) + 4hv{v)] + 6fE{v) + 3hE{v) 



(18) 



This expression is independent of the geometric parameters Oj as expected from Nernst's 
third law of thermodynamics. The direct evaluation of these integrals gives the result 



S 



1 /"^ 1 / 1 

-- / dv-[ coth(t>) — Ky(v) 

4 Jo ^ ' V 

/■°° 1 / 1 

1 + / dv- coth(t;) - Ke{v) 

. n v\ V 



(19) 



where the first term comes from the contribution of the volume terms, and the second from 
the edge terms. Now, if K{v) is taken strictly either as 0, or as 1 or as —1 for any value of 
V, the entropy at zero temperature will not be zero; even worse, it will clearly diverge. If 
K = 0, the integrands would be regular at f = 0, while the edge and volume integrals would 
have a logarithmic divergence. On the other hand, if K equals 1 or -1, the integrand would 
not be regular at f = and the integrals would diverge. Thus, the simplest way of keeping 
a bounded value for the entropy is to choose K = for v < vo and K = 1 for larger values 
of V, with vq to be determined. Accordingly, we define the function 

G(vq) = dv-( cothft;) --)+ [ dv-( cothft;) - - - iV (20) 
Jo v\ v) J^^ v\ V J 

and choose a value of v^^ in such a way that the entropy does take the value zero at T = 0. 
A plot of G'(fo) is shown in Fig. 1; it turns out to be a increasing function of fo, such that 
G[vy) = and G{ve) = -1 for vy = 1.763876988 and ve = 0.64889408 (with ten significant 
figures) . 

Thus choosing Ky^E = 6(f — vv,e), where 6 is the standard step function, we guarantee 
that AF satisfy the boundary condition Eq. (fT^ and that a finite regularized value of the 
thermodynamic variables AF and AS* is obtained. 

At this point, we recall that the spectral density given by Eq. 1^ is not valid for uj = 
0; however, if this equation is taken as it stands for all values of u, it includes a term 
proportional to 6{uj). It is precisely such a term that gives rise to the extra factor 3/4 
appearing in the contributions of the edges to the entropy at zero temperature in Eq. (|19|). 
This unphysical contribution to the entropy is properly removed by the selection of f^;. 

From a physical point of view, it is not possible in general to distinguish between processes 
that differ from each other by the emission of a certain number of low energy photons with 
u ^ 0. In free space, this is true for a continuum, but in a closed box, the values of vv,e 
are the ones that define which photons should be taken as being of low energy; namely, 

Vn>vv^ keT > (21) 

UnTl 

and this condition determines whether the photons contribute to the counting of microstates. 
That is to say, the uniqueness of the vacuum determines which modes are to be considered 
as infrared. 

Given S and F, the expression for the internal energy can be directly calculated as 

E{T) = Ebb{T) + AE{T) , 

9 



where 



and 



A£;(T) = AF + TA5 



010213 

- -^T^Yl [^^d'ivnoo) + aafi-Vono) + asg'ivoon)] . (22) 

n 

Notice that AE{T) does not depend on the value of K. In fact, this expression can be 
obtained directly from the spectral density without any explicit calculation of AS and 
AF. The regular behavior of AE is the reason why previous authors have restricted their 
calculations to this particular thermodynamical variable 

The behavior of the total entropy and the difference between its "free space" and bounded 
space values, is shown in Fig. 2, together with the free and internal energies. Three particular 
configurations were considered for the numerical calculations: the cases of a "pizza box", 
^2 = 03 = lOOai; a cube, ai = 02 = 03; and a wave-guide configuration, 02 = as = ai/10. 
Some common features are worth noticing: (i) the total entropy S and the relative entropy 
AS are zero at T = 0, as it should be by construction; (ii)the total entropy is always positive; 
(iii)the relative entropy may take negative values and, accordingly, heat may be released by 
the system when going from free space to a bounded configuration; (iv) the entropy and free 
energy of the system are discontinuous functions of the temperature; (v) the free and internal 
energies equal the Casimir energy at T = 0; (vi) all relative thermodynamic variables tend 
to zero as T ^ 00. 

The case 02 = 03 = lOOoi is similar to the widely studied parallel plates configuration; it 
must be realized, however, that there is a qualitative difference between the two configura- 
tions: in closed rectangular boxes, all frequencies are discrete and no zero frequency electro- 
magnetic modes are allowed, whereas parallel plates admits a continuous range of frequencies 
and the existence of either TE and TM of zero frequency is not excluded in principle|2fl|. 



The presence of such modes significantly alters the behavior of the free and internal 



energies 



ener 

m 



at high temperatures and introduce terms that increase linearly with temperature 
Our calculations show that for a closed rectangular cavity, such linear terms are absent and 
that, moreover, all the effects associated to the discrete density of modes Ap{uj) tend to 

10 



disappear as the temperature increases. The case of a cube, ai — a2 — as, is particularly 
interesting since its Casimir energy is positive; in this case, the difference between the free 
and bounded internal energy AE is a decreasing function of the temperature. In the case 
Oi = 02 = 03/10, the cavity is similar to a waveguide and discontinuities of F and S persist 
at higher temperatures. 

Once a regularized expression for the thermodynamic potentials has been obtained, other 
relevant physical variables can be directly calculated. For instance, the specific heat Cy for 
a given geometry. 



is clearly a continuous function of T. For a cube, it turns out that the presence of the 
boundaries reduces the specific heat Cy with respect to the case with boundaries at infinity, 
while the opposite effect occurs for the "pizza box", oi <S 02 = 03. As a general feature, 
Cy decreases at very low temperatures with respect to the free space configuration, but 

increases at moderate temperatures. 

The pressure on wall 1 of area 0203 is given by 



02^3 ^daiJ T ' 

with similar expressions for P2 and P3. A straightforward calculation shows that, as ex- 
pected, the equation of state 



is satisfied. For a cube Pi — P2 — P3 — E/3V and the pressure is a continuous functions of 
T . For other configurations, the regularized expression of F leads to discontinuities of the 

pressures at the walls at low temperature. This can be seen in Fig. 3, where a plot is shown 
of the total pressure at walls 1 and 3 for "waveguide" and "pizza box" configurations. It 
is worth recalling that the pressure is the physical variable most accessible to experimental 
verification. 

A common feature of some thermodynamic variables for all configurations is the ap- 
pearance of discontinuities as the temperature varies. These discontinuities are due to the 
different weights that the EM modes acquire as the temperature increases. 




(23) 




E={P, + P2 + P3)V 



(25) 
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III. OTHER APPROACHES 



In this section we compare our results with those reported in the hterature for the same 
or similar systems. For this purpose, it is convenient to use the formula 

1 2x 1 
cothvra; = — > — — (26) 

TVX TT ^ — ^ + /C^ 

k=l 



and write the regularized free energy in the form 

V 



AF = AEn + AFm + 



\Vn\>Vv 



T, , , , , , 

VnO(l>VE VQriO>VE V(iOn>VE 

where AFm is the free energy that follows directly from the Matsubara formalism, see e. g. 
Santos and Tort {2^. It can be obtained in our own formalism using Ap with the integration 
paths so chosen that Ky and Ke are taken as zero for every Vnim- Explicitly: 

2Vh^'^ 1 hai ^ 1 



AF - ^'^^V^V ^ ftai 



U. 



nOO 



/ia2 1 ^fls 1 . . 

Notice that this expression has a logarithmic divergence and, consequently, it precludes any 
practical calculation of physical quantities. The analogue of AFm for a massless scalar field 
was reported in Ref. p, where, as mentioned above in the Introduction, the problem of 
the divergence was stated but no explicit finite expression for the free energy, valid at all 
temperatures, was given. 

From the explicit form of the terms that depend on Kv,e{v) in Eq. ()12|). it can be seen 
that they have a quantum origin: although each term in the summation is independent of 



h, the selection o 
Santos and Tort 



the integers rii that contribute to the summation depends on h. Indeed, 
have shown that the free energy of the EM field calculated using the 
Z-function regularization technique is identical with AFm up to a term T\n(fi / ^/2tt T), 
where /i is a scale factor. Santos and Tort choose /i = \^2n T, whereas in our regularization 
scheme 



lfn|>l'V' ''nOO>'"-E '>'OnO>'"E '!'00n>'!'_B 

12 



is the scale factor. 

Since it is well known that most difficulties with the quantization of the electromagnetic 
field are related to the fact that it is a massless field, the origin of our divergence problem 
can also be clarified by assigning an effective mass to the photon. For a massive field, 
the volumetric contribution to AFm must be changed to p: 

Vmlc^, » A'J(m,cV2R) \{hk/Ty- + ul] 

where K2 is the associated Bessel function. The limit of low temperature, T ^ 0, of this 
expression can be calculated as a continuous integral over variables n^, with the result 

AFM(m^) ~ -T In (l - g-'^-^'/^T j ^ (3^^ 
from where the following contribution to the entropy is obtained: 

which is finite and does tend to zero in the limit of low temperature. Clearly this limit is 
to be understood in the sense that T -C m^c^. Had we taken = from the beginning, 
the entropy would not be finite; thus the origin of the divergence can be traced back to 
the fact that, since the photon is massless, there is no natural energy for the temperature 
to be compared with; this seems to be the origin of the difficulties with the definition of 
entropy in closed rectangular cavities. In principle, one could use Eq. (j30|l to calculate the 
thermodynamic variables and take the limit of zero mass only at the end, but this procedure 
is not useful for practical computations; as far as we know, no expression for the regularized 
quantities in closed form has been obtained in this way. 



IV. SUMMARY AND CONCLUSIONS 



In this paper we have obtained fully regularized expressions for the thermodynamic vari- 
ables associated to the electromagnetic field inside a rectangular box with perfect conducting 
walls. The regularization procedure we have used is based on the condition that the entropy 
be zero at T = 0, which is equivalent to imposing the uniqueness of EM vacuum. In this 
way, we have been able to calculate thermodynamical quantities that could be tested exper- 
imentally. 
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In general, infrared divergences in these type of calculations are due to the possible 
emission and absorption of an indefinite number of soft photons. However, in a closed 
rectangular cavity of maximum size L, the frequencies are discrete and there must be a 
lower bound to the frequency, ujmin = cf^L, such that the limit u; — > is never achieved in 
the system. According to our analysis, this seems to be the origin of the difficulties in the 
infrared limit. Thus special care must be taken in counting the different microstates in the 
evaluation of the free energy: specifically, it is necessary to impose the uniqueness of the 
state without photons of any allowed frequency. Otherwise, the entropy is not well defined 
at zero temperature. 

In the regularization scheme we propose, we have taken the above facts into account 
using a cut-off procedure. According to quantum statistics, the cut-off values vy and ve 
can be interpreted as parameters that define which EM modes are compatible with the 
macrostate of a given configuration, that is, those that must be properly counted at low 
temperatures. In particular, it turns out that the cut-off related to the volume and the 
edges of the configuration are different: vy 7^ ve- 

The introduction of cut-off terms is also related to the fact that all energy fluctuations 
induced by the thermal bath must be larger than the quantum fluctuations. Actually, this 
is the meaning of the important inequality ()21|). On the other hand, if the photons had an 
effective mass ~ h/Lc, no cut-off would be necessary. 

The discontinuities exhibited by the free energy and the entropy, associated with the 
modes that flt into the configuration as the temperature varies, are some of the most striking 
predictions of our calculations. However, they could be due to an excessive idealization of 
the system under study, namely, a perfectly conducting cavity at thermal equilibrium with 
its surroundings. It may happen that for a dielectric closed box such discontinuities are 
softened, but traces of them could be present in experiments. 
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Appendix 



A basic integral to be evaluated is: 

— — / (ia;cos(tia;) In ( 1 — e"^'^ j 
ph Jo ^ ^ 

sm{uu!) 







dcu- 



^ poo poo 

- / dujsm{uuj) coth.{ph(jj/2) — - / dujsm{uuj) 
4 J-oo 2 Jq 



71 



y- (33) 



2ph 

(for K 7^ 0). The value of this integral depends on the integration path used to circumvent 
the pole at a; = 0. For a given u > 0, one finds that X = ±1 if the path circumvents the 
singularity over (under) the real axis, and K — if the principal value of the integral is 
taken. 

The other basic integral is 

--— / dcucusmiuuj) In ( 1 



■-^ / du;u;sm{uio)ln(l-e-^'^) 

ph Jq V / 



I 



dcvuj^ + - / duj-r^ , 34 



^ " Jo 



, C0S(UUJ) 1 , , , , . /n^ 

duj uj—^ — 7 u!cos[uuj) coth[pnLj/2) — 



■00 

00 



duj UJ cos(ww) 



with 



1 

Notice that the value of this last integral is independent of the integration path around the 
origin, a; = 0, because the residue at that point is zero. 
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Figure Captions 



Figure 1. Plot of the monotonic increasing function G{vo). Numerical analysis shows 
that G{vv) = for w = 1.763876988 and G{ve) = -1 for ve = 0.64889408 (with ten 
significant figures). 

Figure 2. Dimensionless thermodynamic potentials S/ks, AS/ks, A/ = naiAF/hc, 
and Au — iraiAE/Hc as functions of the dimensionless variable ^ = TrkBTai/hc. Three 
particular configurations are shown: 02 = 03 = lOOai (pizza box), Ox — 0,2 — '^3 

(cube), and 

/lO (wave guide), in the first, second and third columns respectively. The total 
entropy is given in the first row; notice that it is always positive, which is not the case for 
the difference between the entropy in "free space" and its finite domain value, as seen in 
the second row . The finite domain free and internal energies are shown in the third and 
fourth rows respectively. In all cases, the finite domain thermodynamic potentials have the 
physically expected values as TrkBTai/Hc — > 0, and tend to zero as TrkBTai/hc — > 00. 

Figure 3. Dimensionless pressures on the walls 1 and 3, pi = TrajPi/hc and pa = 
nafP^/hc as function of the dimensionless variable ^ = TikBTai/hc for a "pizza box" config- 
uration (upper panel), and a "waveguide" configuration (lower panel). For comparison the 
pressures corresponding to Fbb are plotted in dotted lines. 
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